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ABSTRACT 

The discovery of multiple stellar populations in globular clusters raises fundamental 
questions concerning the formation and dynamical history of these systems. In a pre- 
vious study aimed at exploring the formation of second-generation (SG) stars from 
the ejecta of first-generation (FG) AGB stars, and the subsequent dynamical evolution 
of the cluster, we showed that SG stars are expected to form in a dense subsystem 
concentrated in the inner regions of the FG cluster. In this paper we explore the impli- 
cations of the structural properties of multiple-population clusters, and in particular 
the presence of the inner SG subsystem, for the disruption of binary stars. We quantify 
the enhancement of the binary disruption rate due to the presence of the central SG 
subsystem for a number of different initial conditions. Our calculations show that SG 
binaries, which are assumed to be more concentrated in the cluster inner regions, are 
disrupted at a substantially larger rate than FG binaries. Assuming a similar initial 
fraction of FG and SG binaries, our dynamical study indicates that the SG popula- 
tion is now expected to contain a significantly smaller binary fraction than the FG 
population. 
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1 INTRODUCTION 

Spectroscopic and photometric observations reveal that 
globular clusters host multiple stellar populations. Spectro- 
scopic studies show star-to-star variations in the abundances 
of light elements, such as Na, O, Al, and Mg, indicating 
that a significant fraction (50-80%) of globular cluster stars 
must have formed out of matter processed through a high- 
temperature CNO cycle in a first generation (hereafter FG) 
of stars (see e.g. Carretta et al. 2009a, 2009b and references 
therein). Photometric studies reveal the presence of multi- 
ple main sequences, subgiant, and red-giant branches in nu- 
merous clusters, supporting the spectroscopic evidence for 
multiple populations (see e.g. Piotto 2009 and references 
therein). Photometric measurements also provide evidence 
for a population of very He-rich stars among the second- 
generation (hereafter SG) stars of some clusters (see e.g. 
Piotto et al. 2007). These findings confirm the predictions 
of previous studies (D Antona et al. 2002; D Antona & Caloi 
2004, DAntona & Caloi 2008) that suggested the existence 
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of a population of stars with a strong He enhancement on 
the basis of the extension of the horizontal branch. 

The origin of the gas out of which SG stars formed is 
still an open question: in addition to the scenario involv- 
ing AGB stars (Cottrell & Da Costa 1981, Ventura et al. 
2001) described in more detail below, possible gas sources 
proposed in the literature also include rapidly rotating mas- 
sive stars and massive binary stars (Decressin et al. 2007, 
De Mink et al. 2009; see also Renzini 2008 and references 
therein for a review). Many fundamental questions concern- 
ing globular cluster star formation and cluster chemical and 
dynamical history are raised by the discovery of multiple 
populations, and are targets of ongoing investigations (see 
e.g. D'Ercole et al. 2008, 2010, Vesperini et al. 2010, Bekki 
2011 and references therein). 

D'Ercole et al. (2008) explored the formation and dy- 
namical evolution of multiple populations in globular cluster 
by means of hydrodynamical and N-body simulations, focus- 
ing on a model in which SG stars form from the ejecta of FG 
AGB stars. Our simulations show that the AGB ejecta form 
a cooling flow and rapidly collect in the innermost regions 
of the cluster, forming a concentrated SG stellar subsystem. 
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In order to form the numbers of SG stars observed today, 
the FG cluster must have been considerably more massive 
than it is now, and the majority of stars in the cluster ini- 
tially belonged to the FG population. The N-body simula- 
tions presented in D'Ercole et al. (2008) show that the early 
expansion triggered by the loss of mass in the form of SNII 
ejecta and primordial gas leads to a strong preferential loss 
of FG stars. 

According to our models, during this early evolution- 
ary phase the cluster evolves from a configuration in which 
FG stars dominate to one in which the numbers of FG and 
SG stars are similar (or even one in which the SG stars are 
now the dominant population), as observed in several Galac- 
tic globular clusters (see e.g. Carretta et al. 2009a, 2009b). 
At the end of this early phase, the SG subsystem, although 
now composed of a number of stars similar to that of the FG 
population, is still concentrated in the cluster inner regions. 
The system thus begins its long-term, relaxation-driven evo- 
lution with a structure characterized by the superposition of 
a compact SG inner cluster and a more extended FG cluster. 
In the initial simulations discussed in D'Ercole et al. (2008) 
we showed that, as the cluster evolves, the two populations 
tend to mix and the characteristic structure imprinted by 
the SG formation process is slowly erased. 

Such a peculiar structure differs significantly from the 
simple King or Plummer models usually adopted as initial 
conditions in studies of globular cluster evolution. It is there- 
fore important to explore the implications of this unusual 
initial state for the long-term evolution of a cluster's struc- 
tural and kinematical properties, as well as its stellar con- 
tent. This paper focuses on the effects of this new class of 
initial structural properties on the survival of binary stars, 
and on the differences in the evolution of the numbers of FG 
and SG binaries. 

Binary stars play a crucial role in cluster evolution, both 
as an energy source supporting the cluster's post-core col- 
lapse dynamics (see e.g. Gao et al. 1991, Goodman & Hut 
1989, Vesperini & Chernoff 1994, Heggie et al. 2006, Trenti 
et al. 2007, Hurley et al. 2007; see also Heggie & Hut 2003 
for a review), and as potential seeds for the formation of a 
variety of exotic objects (e.g. LMXBs, CVs, blue stragglers 
etc.; see e.g. Ivanova et al. 2006, 2008, Ferraro & Lanzoni 
2008 and references therein) . Exploring the evolution of bi- 
naries in the environment resulting from the formation of 
SG stars, and understanding whether or not observed binary 
properties and abundances contain any imprint of a cluster's 
formation, are fundamental problems to be addressed. Here 
we present an initial study of this problem that combines 
the results of a set of N-body simulations of clusters hosting 
multiple populations with analytical expressions for binary 
interaction rates, to estimate the net binary disruption rate 
during cluster evolution. 

The structure of this paper is as follows: in Section 2 we 
describe the N-body simulations and the analytical frame- 
work used for our study; in Section 3 we present our results; 
and in Section 4 we discuss our results and summarize our 
main conclusions. 



2 ANALYTICAL AND NUMERICAL 
FRAMEWORK 

2.1 Binary star disruption 

The outcome of the interactions between a binary and single 
stars in a cluster depends critically on the ratio x = e/ma 2 of 
the binary binding energy, e = Gm 2 /2a to the mean kinetic 
energy of cluster stars, ma 2 , where m is the mass of a single 
star and each of the binary components (assumed here to be 
the same), a is the binary semi- major axis, and a is the 1-D 
stellar velocity dispersion (notice that ma 2 is often also re- 
ferred to as kT in the literature). Substantial analytical and 
numerical effort has been devoted to the study of binary- 
single star interactions and their possible outcomes (see e.g. 
Heggie & Hut 2003 for a review and references therein). 
Binaries having x ^ 1 (so called hard binaries) will, on av- 
erage, become more bound as a result of interactions with 
single stars; for x ^ 1 (soft binaries), binary-single-star in- 
teractions will generally lead to the disruption of the binary 
star. 

The goal of this paper is to explore the extent to which 
the presence of a compact SG subcluster in the inner regions 
of a young globular cluster can increase the disruption rate 
of binaries. We focus here on the process of binary ioniza- 
tion, that is on the disruption of a binary following a single 
interaction with a single star. The disruption rate of a pop- 
ulation of binaries with number density rib, semi- major axis 
a interacting with a population of single stars with number 
density n s can be written 

— -JT = n s ira 2 V th R(x), (1) 
rib at 

where R(x) is an analytical fit obtained by Hut & Bahcall 
(1983) from a series of numerical simulations of binary-single 
star encounters: 

R(x) = (2 ) 

y ' (l + 0.2A/x)[l+exp(x/A)]' { ' 

and Vth = 3(y4/2) 1//2 cr is the dispersion of binary-single-star 
relative velocities. The quantity A is equal to 1 in the case 
of energy equipartition between the binary and the single 
star populations, and equal to 4/3 in the case of velocity 
equipartition. As a cluster evolves toward energy equiparti- 
tion, the appropriate value of A will therefore change from 
4/3 to 1. Hereafter, we adopt A — 1 as our default value. 

In this study we follow the evolution of the overall bi- 
nary disruption rate, and the differences in the disruption 
rates for SG and FG binaries, by combining EqQ]witli the 
time-evolving radial profiles of n s and a obtained from our 
N-body simulations of multiple-population clusters. 

In our calculations we assume that the spatial distri- 
bution of SG (FG) binaries always follows the distribution 
of SG (FG) single stars, n b ,sG = {Nb,sa/N s ,sG)n s ,sG (and 
similarly for FG). We note that, by making this assump- 
tion, we explicitly neglect any possible spatial segregation 
of the binary population and the ensuing increase in bi- 
nary disruption. This neglect will be most significant for 
the initially most compact SG population, and will further 
increase the preferential disruption of SG binaries described 
below. Other processes whose effects are not included here 
are binary-binary interactions and the possible formation 
of binaries by three-body and two-body tidal encounters. 
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Finally, for clusters with a significant fraction of primordial 
binaries, binary heating will delay deep core collapse and the 
binary disruption occurring during that high-density phase. 
(Notice, however, that a larger fraction of binaries delaying 
deep core collapse will increase the number of binary-binary 
interactions, again leading to significant binary disruption). 

In order to quantify the binary disruption rate, we de- 
fine the radial profile of the ionization rate 



I(r, a) = n s (r)%a V t h(r)R(x) 



(3) 



and the volume-integrated disruption rate for SG (FG) bi- 
naries 
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(and similarly for Afg), where a, x, and Vth are defined 
above. 

Finally, we calculate the time integral of Eq[3]to obtain 
the time evolution of the number of SG binaries with a given 
semi-major axis a as 

rt 



N b ,sG(a,t) 
N b ,sG(a,Q) 



— exp 



A SG (a,t)dt 



(5) 



and similarly for Nt,,sG{o>,t). 



2.2 N-body simulations 

As discussed in Section[T] our previous hydrodynamical and 
N-body simulations (D'Ercole et al. 2008) predict that SG 
stars form in a compact subsystem concentrated in the in- 
ner cluster regions (see also Bekki 2011 for additional sim- 
ulations confirming our prediction), and that early cluster 
expansion triggered by the loss of SNII ejecta and primor- 
dial gas is responsible for the loss of a significant fraction of 
FG stars. 

In order to explore the long-term evolution of a 
multiple-population cluster, driven by two-body relaxation, 
we have followed the evolution of four different systems each 
containing equal numbers of SG and FG stars. In all cases, 
the initial state of the FG system is modeled as a single- 
mass King (1966) model with central dimensionless poten- 
tial Wo = 7 and truncation radius R t equal to the cluster 
Jacobi radius. The initial SG subsystem is also modeled as 
a Wo = 7 King model, but one that is concentrated in the 
inner regions of the FG cluster. We have explored the evo- 
lution of the resulting two-population cluster for four val- 
ues of the initial ratio of the FG to SG half-mass radii: 
Rh,FG/Rh,SG = 2.5, 5, 10, 25. We refer to these simula- 
tions as r2p5, r5, rlO, and r25, respectively. The different 
ratios represent clusters with the same initial mass and to- 
tal radius (Rt) but different internal structures, as measured 
by the degree of concentration of the SG subsystem. 

All simulations start with a total number of particles 
TV = 10, 000, and our analysis focuses on the evolution of 
the system until t ~ 40t r h,sG(0), where t r h,SG(0) is the 
initial half-mass relaxation time of the SG subsystem. For 
all initial conditions considered, the system undergoes core 
collapse at t w 17 — 25t r h,sa(0)- More detailed discussion 
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Figure 1. Radial profile of the ionization rate, I(r) (see Eq[3]l 
calculated for a binary with semi-major axis a = 0.05Rt/N for 
the four systems considered, normalized to the central value of 
the SP system. The cyan dashed line shows the profile for the SP 
system; the solid lines show, from top to bottom, the profiles for 
the r25, rlO, r5, and r2p5 systems. 



of the structural evolution of the FG and SG systems will 
be presented in a separate paper (Vesperini et al. 2011, in 
preparation) . 



3 RESULTS 

We consider a common range of binary semi-major axes, 
10~ 2 Rt/N < a < W^Rt/N, for all simulations. To illus- 
trate the effect of the inner SG substructure in our multiple- 
population clusters, in some cases we compare quantities re- 
lated to binary disruption with those calculated for a 'stan- 
dard' Wo = 7 King model with no multiple populations 
and hence no additional substructure due to the presence 
of the SG system (hereafter we will refer to this system 
as the SP -Single Population- system). For the SP system, 
Rh ~ O.IRt and the above range of a corresponds to bina- 
ries that are close to the hard/soft boundary or are hard 
everywhere in the cluster; specifically for the widest binary 
considered, x £ 2.8 in the SP cluster. 

Fig[T]shows the radial profile of the ionization rate I(r) 
for a = 5 x 10 -2 Rt/N in each of the models explored. 
The profiles are calculated using the cluster properties af- 
ter about 5-6 dynamical times, in order to allow the initial 
two-population configuration to come into dynamical equi- 
librium. This figure illustrates how the presence of the inner 
SG system greatly enhances the ionization rate over that in 
the simple SP system. Even for the system (r2p5) with the 
least concentrated SG population, the ionization rate in the 
cluster inner regions is about an order of magnitude larger 
than that in the SP system. 

As the system evolves, the ionization rate reaches a 
maximum at the time of core collapse and decreases again 
in the post-core-collapse phase. However, the ionization rate 
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Figure 2. Time evolution of the volume-integrated ionization 
rate (normalized to the value for the SP system) for SG, Age (a, t) 
(blue lines), and FG binaries, g(a, t) (red lines), for a = 5 X 
10 — 2 Rt/N. Each set (red and blue) of lines refers, from top to 
bottom, to the r25, rlO, r5, and r2p5 systems. 
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Figure 3. Time evolution of the ratio of the SG to the FG 
volume-integrated ionization rate AgQ(a,t)/ApQ(a,t), for a = 
5 X 10~ 2 Rt /N . From top to bottom, the four lines refer to the 
r25, rlO, r5, and r2p5 systems. 



in all multiple-component models is always larger than that 
in the reference SP system. Fig[2] shows the time evolution 
of the volume-integrated ionization rate for SG (Asc(a,i)) 
and FG {A FG (a,t)) binaries, for a = 5 x 10" 2 R t /N. SG bi- 
naries are preferentially located in the cluster inner regions 
and so are disrupted more efficiently than FG binaries. The 
time evolution of AsG(a-,t)/ AFo{a-,t) is plotted in Fig(3] 



Figure 4. Ratio of the number of SG to FG binaries, at t = 
40t r j lj 5G'(0), versus binary semi-major axis, a, for the r25 (red 
dot-dashed line), rlO (black solid line), r5 (blue dashed line), 
and r2p5 (dotted green line). 

As the cluster evolves and the two populations tend 
to mix, AsG(o,,t)/AFo(ci,t) tends to decrease. However, 
the mixing is not complete by the end of the time interval 
spanned by this study (t ~ 4:0t r h,sa(0)), so the SG binary 
disruption rate is always larger than that of FG binaries. 

Fig. Uplots Nb,sG( a )/Nb,Fc( a ), calculated at the end 
of each simulation, as a function of binary semi-major axis. 
In all cases, SG binary disruption is significantly enhanced 
compared to that of the FG binary population. Fig.[5]shows 
the ratios of the final to the initial number of binaries for 
the SG and the FG population, Nb,sG(o) /Nb,sG(a)imt and 
Nb,FG(a>) /Nb,FG(a)init, and illustrates the extent of the bi- 
nary disruption in our multiple population clusters, as well 
as the preferential disruption of SG binaries. 

Finally, in Fig[6] we plot the time evolution of 
Nb,SG(a, t), Nb,FG(a, t), and their ratio, for two different val- 
ues of a for simulation rlO. This figure illustrates the extent 
of the early (pre-core collapse) disruption of binaries due 
to the presence of the high-density SG subsystem. Binary 
disruption further increases during core collapse, and finally 
slows down during the post-core collapse phase. The pref- 
erential disruption of SG binaries continues for the whole 
simulation. 



4 DISCUSSION AND CONCLUSIONS 

The results presented in this paper show that the proper- 
ties of binary stars are significantly affected by the initial 
structure of a cluster hosting multiple populations, and may 
contain important clues to the formation and evolutionary 
history of multiple populations. 

The central result of our study is that significant dif- 
ferences in the numbers of FG and SG binaries is a finger- 
print of the structural properties predicted by our models 
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Figure 5. Numbers of SG and FG binaries (normalized to their 
initial values) at the end of the simulations, versus binary semi- 
major axis, a, for the r25 (red dot-dashed line), rlO (black solid 
line), r5 (blue dashed line), and r2p5 (dotted green line) runs. 
For each set of lines, the upper (lower) curve refers to the FG 
(SG) binary population. 
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Figure 6. Time evolution of the numbers of SG and FG binaries 
(normalized to their initial value) and their ratio, for simulation 
rlO for FG binaries and a = 3 X 10 — 2 Rt/N (upper red solid line), 
SG binaries and a = 3 X 10 -2 Rt/N (lower blue solid line), FG 
binaries and a = 5 X 10 — 2 Rt/N (upper red dot-dashed line), SG 
binaries and a = 5 X 10 — 2 Rt/N (lower blue dot-dashed line). The 
upper black dashed line shows the ratio of the number of SG to 
FG binaries for a = 3 X 10 — 2 Rt/N; the lower black dashed line 
shows the SG to FG binary number ratio for a = 5 X 10 -2 Rt/N. 



of SG star formation and dynamical evolution (D'Ercole et 
al. 2008; further investigation of the long-term evolution of 
multiple population clusters will be presented in Vesperini et 
al. 2011, in preparation). Specifically, we showed in our pre- 
vious studies that SG stars forming from AGB ejecta tend to 
form in a strongly concentrated subsystem in the FG clus- 
ter inner regions (see also Bekki 2011); we have now shown 
in this paper that in a cluster with such a structure, SG 
binaries are preferentially disrupted and that, more gener- 
ally, binary disruption is enhanced compared to a standard 
cluster with similar mass and size but without an inner SG 
subsystem. 

Our calculation is based on analytical calculations com- 
bined with the results of N-body simulations of cluster struc- 
tural and kinematical evolution. Further refinement of the 
calculations presented in this paper will require simulations 
including full treatment of binary stars; simulations includ- 
ing the self-consistent treatment of binaries are very com- 
putationally expensive and are beyond the scope of this ini- 
tial exploratory study. Future studies based on simulations 
with binaries will allow us to incorporate additional effects 
(e.g. binary-binary interactions, binary segregation, binary 
heating) not included in our simple analytical treatment of 
binary disruption. 

The first observational indication of the preferential dis- 
ruption of SG binaries, as found in our analysis, comes from 
a recent study of Ba stars in globular clusters (D'Orazi et al. 
2010). Ba stars are thought to be the result of the accretion 
of matter processed by a thermally pulsing AGB onto the 
secondary component of a binary system. The Ba-rich en- 
velope of the AGB component contaminates the companion 
via wind accretion, possibly followed by stable Roche lobe 
overflow, or by common envelope evolution, depending on 
the mass ratio at the time the donor AGB fills its Roche 
lobe (see Han et al. 1995 for a detailed study of the binary 
channels for the formation of Ba stars; see also Jorissen et 
al. 1998, McClure & Woodsworth 1990 for two observational 
studies of the Ba stars orbital properties). 

In order for this process to occur, the initial binary sep- 
aration must be larger than approximately 1 AU (smaller 
separations would affect the primary component evolution 
before it reaches the AGB phase). D'Orazi et al. find that 
Ba stars belong predominantly to the FG population. This 
result is consistent with a dynamical history in which the 
SG Ba star binary progenitors are disrupted more efficiently 
than those of the FG population. 

We note that the larger binary interaction rate of SG 
binaries can also lead them to harden more rapidly. By tight- 
ening a binary below the minimum separation for the for- 
mation of Ba stars, the more rapid hardening of SG binaries 
may represent an additional channel for the suppression of 
Ba stars in the SG population. Another interesting possible 
consequence of an enhanced binary interaction rate is that 
a binary might be disrupted after the mass transfer episode 
that led to the formation of the Ba star, resulting in a single 
Ba star. 

The sample of Ba stars discussed in D'Orazi et al. is 
still very small, and further observational studies aimed at 
exploring the relative abundance of FG and SG binaries will 
be extremely important to test the predictions of our study 
and, more generally, to shed further light on the interplay 
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between the dynamics of multiple population clusters and 
their binary star content. 

Many observational studies cover only a limited region 
of a cluster, and thus could provide only local information 
on the SG-to-FG binary number ratio. A detailed compari- 
son between observations and theoretical results will entail 
a more detailed accounting of the radial variation of the 
SG and FG binary fractions. The results presented here il- 
lustrate a clear general trend in the evolution of the total 
numbers of SG and FG binaries, but the interplay between 
binary segregation and binary disruption, along with possi- 
ble initial differences in the spatial distributions of FG and 
SG binaries, will presumably result in radial gradients in the 
SG-to-FG binary number ratio. (For example, in the inner- 
most regions, the binary population might be dominated by 
SG binaries even though FG binaries may dominate glob- 
ally.) As discussed above, future simulations including self- 
consistent treatments of binaries will include these effects, 
allowing us to further explore these issues. 



ACKNOWLEDGMENTS 

E.V. and S.M. were supported in part by grants NASA- 
NNX10AD86G, HST-AR-12158.01. F.D and A.D. were sup- 
ported in part by the PRIN-INAF 2009 grant 'Formation 
and Early Evolution of Massive Star Clusters'. 



Heggie, D. C; Trenti, M.; Hut, P., 2006, MNRAS, 368, 677 
Hurley, J. R., Aarseth, S. J., Shara, M. M., 2007, ApJ, 665, 
707 

Hut, P., Bahcall, J.N., 1983, ApJ, 268, 319 

Ivanova, N.; Heinke, C. O.; Rasio, F. A.; Taam, R. E.; Bel- 

czynski, K.; Frcgeau, J., 2006, MNRAS, 372, 1043 
Ivanova, N.; Heinke, C. O.; Rasio, F. A.; Belczynski, K.; 

Fregeau, J. M., 2008, MNRAS, 386, 553 
Jorissen A., Van Eck S., Mayor M., Udry S., 1998, A&A, 

332, 877 
King, I.R., 1966, AJ, 71, 64 

McClure R.D., Woodsworth A.W., 1990, ApJ, 352, 709 
Piotto, G. et al. 2007, ApJ, 661, L53 

Piotto, G., 2009, in Mamajek E., Soderblom D., Wyse, R. 
eds., Proc. IAU Symp. The Ages of Stars, Cambridge Uni- 
versity Press, p. 233 
Renzini, A., 2008, MNRAS, 391, 354 
Trenti, M., Heggie, D. C. Hut, P., 2007, MNRAS, 374, 344 
Ventura, P., D'Antona, F., Mazzitelli, I., Gratton, R., 2001, 

ApJ, 550, L65 
Vesperini, E., Chernoff, D.F., 1994, ApJ, 431, 231 
Vesperini, E., McMillan, S. L. W., D'Antona, F., D'Ercole, 
A., 2010, ApJ, 718, L112 



REFERENCES 

Bekki, K., 2011, MNRAS, 412, 2241 

Carretta, E., et al. 2009a, A&A, 505, 117 

Carretta, E., Bragaglia, A., Gratton, R., & Lucatello, S. 

2009b, A&A, 505, 139 
Cottrell, P. L.; Da Costa, G. S., 1981, ApJ, 245, L79 
D'Antona, F., Caloi, V., Montalban, J., Ventura, P., & 

Gratton, R. 2002, A&A, 395, 69 
D'Antona F., Caloi V., 2004, ApJ, 611, 871 
D'Antona, F., & Caloi, V. 2008, MNRAS, 390, 693 
Decressin, T., Meynet, G., Charbonnel, C, Prantzos, N., 

& Ekstrom, S. 2007, A&A, 464, 1029 
de Mink, S.E., Pols, O.R., Langer, N., Izzard, R.G., 2009, 

A&A, 507, LI 

D'Ercole, A., Vesperini, E., D'Antona, F., McMillan, S. L. 
W., Recchi, S., 2008, MNRAS, 391, 825 

D'Ercole, A., D'Antona F., Ventura P. , Vesperini E., 
McMillan S.L.W., 2010, MNRAS, 407, 854 

D'Orazi, V., Gratton, R., Lucatello, S., Carretta, E., Bra- 
gaglia, A., Marino, A. F., 2010, ApJ, 719, L213 

Ferraro, F.R., Lanzoni, B., 2008, in Dynamical Evolution of 
Dense Stellar Systems, Proceedings of the International 
Astronomical Union, IAU Symposium, Volume 246, p. 
281-290, E. Vesperini, M. Gicrsz, A. Sills Editors, Cam- 
bridge University Press 

Gao B., Goodman J., Cohn H., Murphy B., 1991, ApJ, 370, 
567 

Goodman, J., Hut, P., 1989, Nature, 339, 40 

Han, Z., Eggleton, P., Podsiadlowski, P., Tout, C. A., 1995, 

MNRAS, 277, 1443 
Heggie, D., & Hut, P. 2003, The Gravitational Million-Body 

Problem, Cambridge University Press 



